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ABSTRACT 

It is well established that stellar effective temperatures determined from photometry and spec- 
troscopy yield systematically different results. We describe a new, simple method to correct spectro- 
scopically derived temperatures ( "excitation temperatures" ) of metal-poor stars based on a literature 
sample with —3.3 < [Fe/H] < —2.5. Excitation temparatures were determined from Fel line abun- 
dances in high-resolution optical spectra in the wavelength range of ~3700 to ^7000 A, although 
shorter wavelength ranges, up to 4750 to 6800 A, can also be employed, and compared with pho- 
tometric literature temperatures. Our adjustment scheme increases the temperatures up to several 
hundred degrees for cool red giants, while leaving the near-main-sequence stars mostly unchanged. 
Hence, it brings the excitation temperatures in good agreement with photometrically derived values. 
The modified temperature also influences other stellar parameters, as the Fel- Fell ionization balance 
is simultaneously used to determine the surface gravity, while also forcing no abundance trend on 
the absorption line strengths to obtain the microturbulent velocity. As a result of increasing the 
temperature, the often too low gravities and too high microturbulent velocities in red giants become 
higher and lower, respectively. Our adjustment scheme thus continues to build on the advantage of 
deriving temperatures from spectroscopy alone, independent of reddening, while at the same time 
producing stellar chemical abundances that are more straightforwardly comparable to studies based 
on photometrically derived temperatures. Hence, our method may prove beneficial for comparing 
different studies in the literature as well as the many high-resolution stellar spectroscopic surveys that 
are or will be carried out in the next few years. 
Subject headings: stars: fundamental parameters — stars: abundances — stars: Population II 
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1. INTRODUCTION 

Determining the atmospheric parameters of a star is 
fundamental to characterizing and understanding its na- 
ture and evolutionary state. These parameters are the 
effective temperature Tcff, surface gravity logg, metal- 
licity [M/H] and for ID abundance analyses also the mi- 
croturbulent \rnic- Particularly, the temperature can be 
determined with several methods, such as from photo- 
metric colors, flux calibrated low-resolution spectra, the 
shape of the Balmer lines in a high- resolution stellar spec- 
trum and through forcing no trend of Fel absorption 
line abundances with the excitation potential of the lines 
("excitation balance" ), also from high-resolution spectra. 
The most common technique is to employ broadband 
colors, with many color-temperatur e calibrations avail- 
able for main-sequence stars (e.g., lAlonso et al.l 119961: 
Casagrande et al.l I2010D and giants (e.g.. lAlonso et all 



1999H based on the infrared flux method (IRFM). Pre- 
cise photometry, reliable reddening values and metallicity 
information are necessary when using this method. Us- 
ing fiux-calibrated low-resolution spectra to measure the 
strength of the Balmer jump and comparing the spectral 
shape to grid s of synthetic flux spectra is less common 
(|Bessellll2007l) . 

^ This paper includes data gathered with the 6.5 meter Mag- 
ellan Telescopes located at Las Cainpanas Observatory, Chile. 
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In the absence of photometry or stellar libraries, or to 
test systematic temperature uncertainties, temperatures 
can also be obtained from high-resolution spectroscopy 
alone. Fitting the shape of the Balmer lines with syn- 
thetic spectra of known stellar parameters yields tem- 
peratures, especially for the warmer main-sequence stars 
and subgiants (e.g., iBarklem et al.|[2002[ ). Finally, an- 
other commonly used method is to employ Fel lines 
throughout the stellar spectrum and force their abun- 
dances to show no trend with the excitation potential 
of each line. Compared to using photometry to derive 
temperatures, no reddening information is needed for 
this technique or the Balmer line fltting method. But 
in turn, high-resolution spectra with resolving power of 
R > 15,000 are required, preferably with large wave- 
length coverage throughout the optical range. A list of 
Fe absorption lines covering a large range of excitation 
potentials is also needed. 

The different methods are known to produce tempera- 
tures with systematic offsets from each other. The "ex- 
citation temperatures" are thereby known to yield lower 
effective temperatures than photometry-based color- 
temperature calibrations, often by a few hundred de- 
grees. S tars on t he upper red giant branch are most 
affected. iJohnsonI (j2002) found their 23 metal-poor gi- 
ants to be 100 to 150 K cooler when using the excitation 
method over several color-temperat ure calibration s. Sim- 
ilar results have been reported by ICavrel et al.l (2004), 
lAoki et al.' ('2007) ' LaTetall (|200g ). iFrebel et all (pOlQl . 
HoUek et al. (2011), and many others. The approaches 
and dealings with this issue have been varied, though. 
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iCavrel et al.l (|2004[ ) reached good agreement between 
their photometrically derived temperatures and excita- 
tion temperatures after they excluded strong lines with 
excitation potential of EP = OeV. These often appear 
to have higher abundances than lines with higher exci- 
tation potential. As a consequence, forcing no trend in 
line abundances with excitation potential can yield ar- 
tificially cooler temperatures. Whil e a similar approach 
(cutting EP < 0.2 eV) worked for iCohen et~a l. (2008), 
it d id not yield t he same outcome for iLai et al.l ([20081) 
and iHollek et al.l (J20li). With or without cutting low- 
EP lines, the latter two studies found significant trends 
of abundance with excitation potential to remain. Con- 
sequently, diff'erent authors have adopted either spectro- 
scopic temperatures, or the photometric ones, or tried 
to find a middle ground between the results of the two 
techniques. 

The situation is further complicated by the fact that 
different model atmosphere codes have been used over 
time, producing different results leading to and possibly 
increasing the differences between spectroscopic and pho- 
tometric temperatures. For example, using a model code 
that t reated scattering as true absorption ICavrel et al.l 
(|2004t ) found larger discrepancies than when using codes 
that treated scatteri ng as Rayleig h scatt ering. 

However, recently, IHollek et al.l ()20lH ) carried out an 
abundance analysis of 16 metal-poor giants using differ- 
ent versions of the MOOG code (Sncdcn 1973). They 
found that if scattering is properly treated, the tem- 
peratures were generally even lower than in an older 
MOOG version that treated scattering as true absorp- 
tion - the opposite of what was found by ICavrel et al.l 
()2004[) (although they make no quantitative statement 
about the magnitude of the effect and how it compares 
to their proper-scattering treatment when not excluding 
lines with £;P< 1.2eV). 

This study aims at empirically resolving the issue of 
producing lower spectroscopic excitation temperatures 
when employing the latest version of the widely-used 
MOOG code that now properly treats the scattering 
(jSobeck et al.l 120111 ). We provide simple corrections to 
the slope of the Fel abundances as a function of exci- 
tation potential for metal-poor stars that bring the ex- 
citation temperatures in rough agreement with photo- 
metrically derived values. This way, we can use "the 
best of both worlds" : not depending on photometry, 
reddening and color-temperature relations, while obtain- 
ing temperatures and hence chemical abundances that 
can be easily and relatively fairly compared to results 
based on photometric temperatures. Our pragmatic ap- 
proach does not resolve the underlying processes that 
likely cause the temperature discrepancies which seem 
to be manifold (scattering treatments in codes, line for- 
mation of low and high-E'P lines, wavelength coverage 
of lines and data quality, effects due to non-local ther- 
modynamic equilibrium, etc.) but at least it provides a 
means to facilitate a more homogeneous analysis of stars 
that have temperatures derived from different methods. 
Much work has been invested in improving photomet- 
ric tem perature scales. For exa mple, using the IRFM 
method, iCasagrande et al.l ()2010D showed this photomet- 
ric temperature determination method to be hardly sen- 
sitive to theoretical model parameters, especially the sur- 
face gravity and metallicity (particularly for extremely 



metal-poor stars) . The uncertainty with respect to their 
zero-point (well-calibrated, against interferometric angu- 
lar diameter measurements and HST spectroscopy) using 
solar twins is 15 K, whereas their final IRFM effective 
temperatures have typical uncertainties ranging from 60 
to 90 K for metal-poor stars and slightly lower values for 
more metal-rich stars. The independence of metallicity 
when deriving temperatures this way makes the IRFM 
method, in principle, the preferred one for low-metallicity 
stars. The only major caveat would be that the IRFM 
method sensitively depends on very accurate reddening 
values, as an 0.01 mag change in E{B — V) leads to ±50 K 
change in the final temperature. 

Given these facts, it appears timely to bring the some- 
what plagued excitation temperatures onto a comparable 
scale. Another reason for such a "re-calibration" is that 
the stellar surface gravities, usually derived through the 
Fel- Fell ionization equilibrium, depend on the derived 
temperatures. With lower temperatures, lower gravi- 
ties are derived. Hence, in many cases of cool giants, 
much too low gravities have been obtained (compared 
with isochrones and parallax-derived values), to a point 
that they become unphysical. This effect is particularly 
pronounced in metal-poor stars with lower S/N spectra 
w hich leave few Fell to be measured (e.g., see Table 5 
in IHollek et al.l (|2011f) ) . In the same way, if only pri- 
marily strong lines are measurable in the spectrum, the 
microturbulent velocity becomes very large, often past 
2.5kms~^. Especially, metal-poor stars in dwarf galax- 
ies have been affected by this, since they are all bright 
giants and the data quality is often rather poor given 



their faint apparent magnitudes (JMcW illiam et al.j 
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In summary, adjusting the spectroscopic temperatures 
leads to improved stellar parameters, at least in the 
metal-poor regime and particularly for metal-poor gi- 
ants with [Fe/H] < —2.0. Broadly adopting this scheme 
may prove beneficial for the many high-resolution spec- 
troscopic surveys (Gaia-ESO, GALAH) that will soon 
provide large amounts of stellar spectra that have to be 
analysed homogeneously and consistently. 



2. "CALIBRATION" 



SAMPLE AND ABUNDANCE 
ANALYSIS 



We chose seven metal-poor stars with metallicities of 
—3.3 < [Fe/H] < —2.5 and spanning the largest pos- 
sible temperature range from 4600 to 6500 K covering 
warm main-sequence to cool giants to investigate sys- 
tematic differences between the excitation temperatures 
and photometric temperatures. 

All high S/N, high-resoluti on spectra were obtai ned 
with the MIKE spectrograph (jBernstein et al.ll2003f ) on 
the Magellan-Clay telescope at Las Campanas Observa- 
tory between 2009 and 2011. An 0'.'7 slit width was used 
resulting in a resolving power of i? ^ 35, 000 in the blue 
and ^ 28, 000 in the red spectral range. MIKE spectra 
cover nearly the full optical wavelength range of 3350- 
9100 A, ahhough the S/N is low below 3700 A. The re- 
duced spectra were normalized, the orders merged and 
then radial velocity corrected. The S/N ratio of the fi- 
nal spectra ranges from 70 to 300 at '^ 4100 A, and from 
70 to 500 at - 6700 A per pixel. 

Equivalent widths were measured by fitting Gaussian 
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profiles to the absorption l ines. For a met a l line list we 
used the list described in iRoederer et al.l ()2010D . Be- 
tween ~50 and ~230 Fel and 2 and 26 Fell lines were 
measured in each sample star. See Table [T] for the lines 
used and their measured equivalent widths, although we 
only list lines with reduced equivalent widths of > —4.5. 
In Figure [T] we furthermore show a comparison of the 
~ 200 equivalent width measure ments of HD l 22563 com- 
mon to our lines and those of lAoki et all ()2007[ ) and 
iCavrel et al.l (|2004| ). Not all our lines were measured by 
these two studies, but those that are in common agree 
very well with each other. The agreements of our mea- 
surements are 0.20±0.16mA and 0.25±0.28mA, respec- 
tively, so no significant offsets exist between the three 
data sets. 



-< 100 






20 40 60 BO 100 

EWhd12Z563 this study [mA] 

Fig. 1. — Comparison of equiv alent width measurements for 
HD122563 of this studyw ith those of lAoki et ahl l|2)07} (black open 
circles) and ICavrel et al.l l|2004l ) (red filled circles). The agreement 
is excellent. 



For our ID LTE abundance analysi s, we u sed model 
atmospheres from ICastelli fc Kuruc3 ()2004D and the 
MOOG analysis code of lSnedenI (|1973[ ). albeit the latest 
version that appropriately treats scattering as Rayleigh 
scattering an d not as true absor ption as done in previ- 
ous versions (jSobeck et al.|[20Tll ). The metal line abun- 
dances are presented in Table [T] Our study about stellar 
parameter determinations is based on equivalent width 
measurements of Fe lines between 3700 and 7000 A. 

3. ADJUSTING THE "EXCITATION TEMPERATURE" 
SCALE 

We used our sample of well-studied stars to produce 
a straight-forward method of adjusting the excitation 
temperature scale in order to arrive at spectroscopic 
temperatures that agree with photometric ones. After 
measuring the equivalent widths we first determined the 
stellar temperatures by forcing no trend of the abun- 
dances of individual Fe I lines as a function of their ex- 
citation potential. In the process we varied the micro- 
turbulent velocity to produce no trend of the line abun- 
dances with reduced equivalent widths. Simultaneously, 
to determine the surface gravity, the Fell abundance 
was matched to the Fel abundance. Finally, another 



input parameter of the model atmosphere is the metal- 
licity [m/H] of the star. We use a general prescription of 
[m/H] = [Fe/H] -1-0.25 to account for the overabundances 
of a-elements of [a/Fe] ^ 0.3 as well as small overabun- 
dances of carbon and/or oxygen ([C,0/Fe] ^ 0.3) which 
are all typical for metal-poor halo stars. 

We then carried out an extensive literature search 
to assess the range of photometric temperatures de- 
termined for each star. Table [2] shows a represen- 
tative range of values found in the literature, includ- 
ing the colors and the color-temperature relations used 
for their calculation. For comparison, we also include 
Teff values calculated using the IRFM method in recent 
works whose color- temperature relations are commonly 
used (e.g., Alonso et al. 1996, 1999, Casagrande et al. 
2010). For two stars, CD-24 17504 and G64-12, we 
found relatively few photometrically derived values, and 
most where done using older temperature calibrations. 
We thus calculated photometric temperatures ourselves, 
based on liter ature photometry and r eddening values 
and using the Ca sagrande et all ()2010[ ) scale. Finally, 
for HE 1523—0901, only one photometric temperature 
is available (Frebel et al. 2007a) but at least it was an 
average based on six different colors. 

As can been seen, temperature ranges of several hun- 
dred degrees were found for some cases, making it dif- 
ficult to assess which individual temperatures would be 
"the best" . Hence, we resorted to choosing "common 
sense" values, as our original aim was to re-calibrate our 
excitation temperature scale to a level common with that 
of photometric temperatures found using the more recent 
color-temperature relations. Our initial excitation tem- 
peratures and associated adopted literature values (based 
on photometry) are listed in Table |31 We stress that 
the photometric temperatures we have adopted are not 
strictly averages of the values listed in Table [51 rather 
they are appropriate near-median quantities from the 
range of values found by multiple authors using differ- 
ent color-temperature relations and photometry. 

Figure[2] (upper panel) then shows the differences in ef- 
fective temperatures between our initial excitation tem- 
peratures and our adopted literature values as a function 
of the initial temperatures. We find a strong, linear re- 
lation, with the largest differences at the cooler temper- 
atures of upper red giant branch stars. An (unweighted) 
least-square regression yields the following relation: 

Toff, corrected = Teff, initial — 0.1 X Toff, initial + 670 (1) 

This temperature dependent behavior is completely in 
line with the vast majority of previous studies that re- 
ported differences between spectroscopic and photomet- 
ric temperatures of several hundred degrees, as most an- 
alyzed star were giants. However, main-sequence stars 
are hardly affected, as we find essentially the same tem- 
peratures. We also overplot nominal error bars of 134 K 
which refiect typical uncertainties in temperature deter- 
minations (60 K for photometric and 120 K for spectro- 
scopic values). 

Since we adopted common sense representative photo- 
metric temperatures from the literature that are based 
on a variety of colors and color-temperature relations, we 
also attempted to calculate temperatures ourselves, for 
an exemplary color, V — K, and using thc iAlonso et al.l 
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TABLE 1 
Equivalent Widths and Metal Line Abundances of Stars in the Calibration Sample 



A Species 

[A] 



X 
[eV] 



logs/ 

[dex] 



HD122563 
EW loge(X) 

[mA] [dex] 



HE 1523-0901 
EW loge(X) 
[mA] [dex] 



BD-18 5550 
EW loge(X) 
[mA] [dex] 



CS22892-052 
EW loge(X) 
[mA] [dex] 



HD140283 
EW loge(X) 
[mA] [dex] 



CD-24 17504 
EW loge(X; 
[mA] [dex] 



5889.950 


Nal 


0.00 


0.108 


196.64 


4.00 


175.91 


3.54 


141.05 


3.52 


153.99 


3.72 


112.99 


3.84 


29.91 


2.7( 


5895.924 


Nal 


0.00 


-0.194 


166.54 


3.89 


155.12 


3.51 


117.04 


3.37 


126.36 


3.53 


89.59 


3.70 


19.75 


2.7f 


3829.355 


Mgl 


2.71 


-0.208 


















132.91 


5.17 


85.21 


4.65 


3986.753 


Mgl 


4.35 


-1.030 


33.86 


5.22 






19.58 


4.99 






15.29 


5.26 







(|1999( ) and lAlonso et aH (|1996( ) calibrations. We chose 
V — K since it is least affected by reddening. The red- 
dening in the dire ction of each object w as obtained from 
the dust maps of iSchlegel et al.l (|1998[ ) maps. However, 
it is known that these maps overpredict the reddening, 
in some cases quite significantly. Accordingly, different 
reddening correction were invokedj_Jbllowing the proce- 
dures of Bonifacio et al. (20Q3) and|Mclcndez & Asplunil 
(HOOS). Temperatures were then calculated based on 
the two reddening estimates to investigate the impact 
of these reddening uncertainties. In a few cases, we 
adopted yet a different value. For HP 140283 , we als o 
adopt E{B - V^) = 0.01, following iR^an et alJ (|1996bf ). 
rather than E{B — V) = 0.16 from the dust maps. For 
HD122563, we also adopt E{B - V) = 0.0, following 
IHondaet all ([200l . 

We also find a strong relation between the difference 
in temperatures as a function of the initial spectro- 
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Fig. 2. — Differences in effective temperatures between our initial 
excitation temperatures and adopted literature photometry-based 
values (top panel) and those calculated from V — K (bottom panel). 
We also show the slope from the top panel in the bottom panel (red 
line) The differences are most pronounced at the cooler tempera- 
tures of upper red giant branch stars, whereas main-sequence stars 
are hardly affected. 



scopic temperatures, as can be seen in Figure [2] (bottom 
panel). The slope is even steeper because V—K produces 
temperatures for giants that are warmer than what we 
adopted. Within the error bars, however, the slope of 
our main temperature comparison is still in rough agree- 
ment with what we determined based on V — K colors. 
We note that this comparison only contains five of the 
seven sample stars. We had to exclude two stars (BD— 18 
5550 and G64— 12) because their V—K temperature were 
rather different (more than 200 K in the case of BD— 18 
5550) from what one would expect the temperatures of 
these stars to be, regardless of any employed method. 

In the end we adopt the slope of our initial com- 
parison because our overall aim is to adjust the spec- 
troscopic temperatures to average photometrical values. 
Even with photometry and color-temperature relations 
available, the temperatures derived from different colors 
often vary 100-200 K, as can be seen in Table [51 Find- 
ing a different slope based on the V — K temperatures 
appears to reflect this dispersion. 

We then turned the temperature differences from Fig- 
ure [2] into values for the slope of the abundance trends 
to adopt in the Fel line abundance vs. excitation poten- 
tial diagnostics. This is shown in Figure |31 In the top 
panel we show how the values of the slopes of the Fel 
line abundances after only correcting for the temperature 
and before adjusting the microturbulent velocity and the 
gravity. To guide the eye, we fit the data points with a 
second-order polynomial, which is also shown in the Fig- 
ure. We estimate typical uncertainties in the slope to be 
0.005 dex/eV or less. 

For near main-sequence stars, the spectroscopic 
method delivers nearly the same temperatures as pho- 
tometry, hence there is no significant slope of line abun- 
dances with excitation potential (~ O.OOdex/eV). For the 
giants, however, applying significant temperature correc- 
tions according to Eq. (1) produces an initial slope of 
up to —0.09 dex/ eV. We then carried out the remain- 
ing parts of the analysis: re-determining the microturbu- 
lent velocity and the surface gravity while ignoring the 
slope in the temperature diagnostic plot of line abun- 
dance vs. excitation potential. While the gravity was 
adjusted by up to ^ 0.7 dex for the most extreme cases, 
the microturbulent velocity decreased by less (~0.1 or 
less in most cases). The model metallicity [m/H] was 
usually increased by about 0.1 to 0.2 dex following the 
increase in the Fe abundance. 

Through the process of adjusting the remaining stel- 
lar parameters (gravity, microturbulent velocity, metal- 
licity), we actually found that the slope of abundance 
with excitation potential was reduced. Figure [3] (bot- 
tom panel) shows the final slopes in the temperature 
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TABLE 2 

Effective Temperatures Collected from Literature 



Trf 



Colors 



Method 



Reference 



Comments 



HD122563 



4670 


28^ 
302 
61 

44I 

75I 
171 


V,J,K 

B,V,R,I 

IRFM 

B,V,R,J,K 

B,V,K 

B,V,K,b,y 

B,V,R,J,K 


COH78, FRO 83 


4650 


BG78, BG89 


4572 


their calibration 


4617 


AL099 


4546 


AL099 


4616 


AL099 


4843 


RM05 



Peterson ct al. (1990) 
Ryan et al. (1996b) 
Alonso e t al. (1999) 
Cavrel et al. (2004) 
Honda et al. (2004) 




phot, from NIC78, CAR83 

their phot., BN84, GS91 

as quoted in PASTEL database 

phot, from AL098, BEE07, EPC99; adopt 4600 

phot, from Simbad; adopt 4570 

phot, from Simbad, HM98 

phot, from CAY04, HON04 



HE 1523-0901 



B,V,R,I,J,H,K AL099 



Frebel et al. (2007a) ' 



4630 



40 



phot, from BEE07 



BD-18 5550 



4750 
4789 
4668 
4700 
4558 



0^ 

271 

129^ 



V,J,K 

B,V,R 

IRFM 

B,V,R,J,K 

B,V,R,J,K 



COH78, FR083 
BG78, BG89 
their calibration 
AL099 
RM05 



Peterson et al. (1990 ) 
McWilliam et al. (19951 
Alonso et al. (19997 ^ 
Cavrel et al (2004) 
Yong et al. (20137 



phot, from CAR78, CAR83 

their phot. 

as quoted in PASTEL database 

phot, from AL098, BEE07, EPC99; adopt 4750 

phot, from CAY04 



CS 22892-052 



4763 
4850 
4860 
4740 
4825 



541 
83I 



B,V,R 

B,V 

B,V,R,I,J,K 

B,V,K 

B,V,R,J,K 



BG78, BG89 
sec RYA96b 
AL099 
AL099 
RM05 



McWilliam et al. (1995) 
Norris et al. (1997) 
Cavrel et al. (2004) 
Honda et al (2004) 
Yong et al. (2013) 



their phot. 

phot, from RYA89, BPS92, 

phot, from AL098, BEE07, 

adopt 4790 

phot, from CAY04, HON04 



NOR97 

EPC99; adopt 4850 



HD140283 



Alonso ct al. (1996) 
Bikmacv ct al. (1996) 
Rya n et al. (1996b) 
"nohen et al. (2002) 
Honda et al. (2004) 
Jonsell et al. (2005) 
^asagrandee^^aL ^2010) 
Yonget al. (2013) 



5691 
5640 
5750 
5750 
5609 
5751 
5777 
5711 



50^ 

53 

34I 

50^ 

55 



IRFM 

V,R,I,K 

B,V,R,I 

V,J,K 

B,V,K 

uvby-/3 

IRFM 



their calibration 

their calibration 

B086, BK92 

HOUOO 

AL096 

AL096 

their calibration 

CASIO 



as quoted in HOS09 

as quoted in MAS99 

their phot. 

their phot. 

phot, from Simbad; adopt 5630 

phot, from OLS83, SN88 

phot, from COH02, HON04 



CD-24 17504 



6060 

6373 
6237 
6455 
6236 
6290 



54-^ 



791,3 
62 



B,V,R,I,b,y 
IRFM 

b,y 
IRFM 



B086, BK92, MAG87, 

VB85 

their calibration 

CAR83, KIN93 

their calibration 

CASIO 

CASIO 




^^ai^_e^^L 

Alonso et al. (1996) 
Primas ct al. (200Qb) 
Cas agrande ct al. (2010) 
Yong et al. (2013) 
This Study 



phot, from RYA89 

as quoted in PRIOOb 
phot, from RYA99 

phot, from RYA91 
phot, from RYA91 



G64-12 



BUS87, B086, KUR89, 

MAG87, VB85 

B086, BK92, MAG87, 

VB85 

their calibration 

AL096 

IRFM 

CASIO 



Rvan et al. (1991) 
Rvan et al. (1996a) 



phot, from RYA89, SN88, CAR83, this work 

phot, from CAR78, CAR89, CS92, SN89, SCH93 

as quoted in PRIOOa 

CASIO 

phot, from AOK06 colors 



6325 
6290 



44 



B,V,R,I,b,y 
2 B,V,R,I,b,y 



6470 90 IRFM 

6430 75I B,V,R,I,K 

6464 61 

6450 35I B,V,R,I,K 
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Aoki et a l. (2006) 
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HOUOO -[tSudashclt et al.l f20 00|); HO N04 -[Honda et al. 120041); HOS09 
ct al. I 1991); MAG87 - Magain ^l9871): MAS99 - Mashon 

et al. 12010); PRIOOa - Prim as et al J l2QQ0a ) ; PRIOOb _ _,, 

RYA96b - Ryan et al I 1996b ); RYA99 - |Uyan ct al.. 1:1999') ; SN88 -Schuster Jk NissenI jlQ^^ ); SN89 
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the average temperatures (column 1) calculated from the different colors (colu 
lent of errors in photometry /reddening. 
. two different color-T„ff relation. 
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TABLE 3 

Stellar Parameters of the Sample 



Star TcB,!!idop. T'eff,V-K 

[K] [K] 



'-^off 

[K] 



initial determination 



[dcx] 



W7W 
[dcx] 



[knis 



after temperature correction 

'-'off 



[K] 



iog(Sj 
[dcx] 



WJW 
[dex] 



[kms 



HD122563 


4650 


4644 


4380 


0.10 


-2.93 


2.65 


4612 


0.85 


-2.79 


2.30 


HE 1523-0901 


4630 


4688 


4380 


0.05 


-2.97 


2.85 


4612 


0.80 


-2.85 


2.70 


BD -18 5550 


4720 




4570 


0.80 


-3.20 


2.05 


4783 


1.25 


-3.02 


2.00 


CS22892-052 


4800 


4903 


4620 


0.85 


-3.24 


2.20 


4828 


1.35 


-3.08 


2.15 


HD140283 


5700 


5637 


5550 


2.95 


-2.74 


1.50 


5665 


3.20 


-2.64 


1.45 


CD -24 17504 


6280 


6276 


6210 


3.60 


-3.26 


1.35 


6259 


3.65 


-3.23 


1.40 


G64-12 


6420 




6430 


4.35 


-3.30 


1.45 


6457 


4.40 


-3.28 


1.50 
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diagnostic plots. The most extreme one is now only 
— 0.07dex/eV which, in the end, only corresponds to 
a temperature difference of 150 K compared to a slope 
of 0.000 dex/eV. However, we caution that if one would 
change the temperature again to produce no trend, the 
gravity and microturbulent velocity would change again, 
and one would once again iterate towards the lower 
temperature/lower gravity/higher microturbulent veloc- 
ity solution. Regardless, 150 K is of the same order as 
what uncertainties in spectroscopic effective tempera- 
tures generally are. It is thus reassuring that despite 
of applying our temperature corrections, we are produc- 
ing temperatures that nominally are no more than 150 K 
away from a trend with no slope. Table |3] shows the ini- 
tial and corrected final sets of stellar parameters of all 
sample stars. 

In Figure IH we show the Fe I line abundances as a 
function of excitation potential for all sample stars. Plot- 
ted are the line abundances showing no trend of excita- 
tion potential as initially derived and the final line abun- 
dances after temperature correction and subsequent ad- 
justment of v„iic and logg. This illustrates how the slope 
changes as a result of our analysis, which is also quanti- 
fied in the two panels in Figured 

We then checked how the initial and corrected sets of 
stellar parameters compare with an isochrone. We used 
12Gyr isochrones with [Fe/H] = -3. 0. -2.5 and -2 
and a-enhancement of [a/Fe] = 0.3 (jKim et al.l 120021 ). 
Figure [5] shows the results. The initial, too cool temper- 



0.02 

-0.02 
-0.04 
-0.06 
-0.08 

-0.1 



0.02 

-0.02 
-0.04 
-0.06 
-0.08 

-0.1 




After temperature correction 



5000 5500 6000 6500 

^ eff .spectroscopic 




After final parameter adjiistments 



4500 5000 5500 6000 6500 

^ eff .spectroscopic 

Fig. 3. — The slopes of the Fel line abundances as a function 
of their excitation potentials after only adjusting the temperature 
according to our method (top panel) and after also adjusting the 
surface gravity and microturbulent velocity (bottom panel) to ar- 
rive at the final stellar parameters. The maximum slope in the 
bottom panel nominally corresponds to ~ 150 K compared with a 
slope of 0.000 dex/eV. 



atures pushed the giants off the isochrones at < 4500 K. 
The corrected values agree very well with the isochrones. 
The only small exception is the subgiant HD140283, 
which is ^0.3 dex off the subgiant branch or 250 K from 
the base of the giant branch. However, the initial stellar 
parameters made this star more of a warm giant rather 
than a subgiant. Hence, the corrected stellar parameters 
provide better solution for this star. 

Finally, we checked how the microturbulent velocities 
as a function of the surface gravity compare with those 
of other studies. In Fi g ure El we co mpare our valu e s with 
those of iCavrel et all (|2004f l and rBarklem et al.l (|2005l ) 
who used photometrically determined temperatures for 
their studies. There exists some scatter in the micro- 
turbulent velocities for a given gravity within each study 
and perhaps also some small systematic difference. While 
the Cayrel et al. values have been derived with the same 
methodology as used in this study, the Barklem et al. 
values are derived from synthesizing spectral lines by si- 
multaneously solving for all four stellar parameters (al- 
though based on app ropriate initial guesses). Table 3 of 
I Barklem et al.l (2005") suggests that with this technique 
slightly higher surface gravities and higher microturbu- 
lent velocities (0.1 to 0.3kms~^) are derived, compared 
to other literature studies. Assuming typical uncertain- 
ties of 0.3kms~^ our microturbulent velocities are in 
good agreement with those from the literature. 

Based on the properties of our chosen sample of stars, 
this method has been accurate for metal-poor subgiants 
and giants with —3.3 < [Fe/H] < —2.5 that have high- 
resolution spectra with a large wavelength coverage (at 
least 4750 to 6800 A). 

We carried out further tests to asses the validity of our 
adjustment scheme at [Fe/H] < —3.5. As pa rt of this we 
analyzed CD-38 245 ([Fe/H] ^ -4.2; e.g. . lYong et al. 
2013). HE 1300-^0157 ([Fe/H] = -3.7; i Frebel et al. 
2007b), HE 0557-4840 ([Fe/H] = -4.8; 'Norri s et al. 
2007), HE 0107-5240 ([Fe/H] = -5.4; Christlie b et al. 
2002), and HE 1327-2326 ([Fe/H] = -5.6: IFrebel et al. 
2005). In all cases, we found that the excitation temper- 
atures yield much lower values ( AToff of 500 to 1400 K) 
for these ultra metal-poor stars than what was found 
for our sample of seven extremely metal-poor stars (the 
maximal AToff is 270 K, see Figure [2]). Hence, we spec- 
ulate that at [Fe/H] < —3.5 the excitation method to 
determine effective temperatures may not be sufficient 
for the task, perhaps because of (metallicity dependent) 
effects, such as NLTE effects or the accuracy of model 
atmospheres (e.g. plane-parallel models for cooler, ex- 
tended giants), or simply because at those metallicities, 
only relatively few Fe lines are available from the spectra. 
While we used all known stars with [Fe/H] < —4.0 that 
had published equivalent widths, the number of stars re- 
mains small, and our individual ATeff values should per- 
haps be regarded as somewhat anecdotal. In summary, 
we suggest that our calibration appears to be breaking 
down at [Fe/H] 4.0. 

We also tested if any of our results would change if 
fewer blue lines would be available due to the use of spec- 
trographs other than MIKE on Magellan or Xshooter on 
VLT that yield (nearly) full optical coverage. Mimicking 
a Subaru/HDS two-arm setting (4000 to 6800 A) and the 
VLT/UVES 580nm setting (4750 too 6800 A), we checked 
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Fig. 4. — Fel line abundances as a function of excitation potential x iov all sample stars. Plotted are the line abundances showing no 
trend of excitation potential as initially derived (blue squares), and the final line abundance trends following the temperature adjustment 
and after Vmic^ logg and metallicity have been adjusted (red triangles). Shown also are the corresponding initial (in brackets) and final 
temperature. 
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- Y2 isochrones, t = 12 Gyr 
[Fe/H] = -3.0, -2.S, -2.0 




6500 



6000 



5500 



5000 



4500 



4000 



^ eff, spectroscopic, final 



Fig. 5. — Our initial (small black open circles) and corrected 
stellar parameters (larger open blue circles) plotted over 12 Gyr 
isochrones with [Fe/H] = -3.0, [Fe/H] = -2.5 and [Fc/H] = = -2 , 
all having an Q-enhancement of [o/Fe] = 0.3 I IKim et al.ll2003) . 
The agreement of the stellar parameters with the isochrones fol- 
lowing the temperature adjustment is very good. The cooler giants 
do not extend past the upper red giant branch of the isochrones 
anymore. 
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Fig. 6. — Final microturbulent velocities as function of the sur- 
face gravity for the sample stars (black filled circles ) as well as the 
ICavrel et al.l l|2004l ) sample (small red circles) and IBarklem et al.l 
1)20051 ) sample (without the horizontal branch stars; small blue cir- 
cles). 



whether using reduced wavelength ranges would change 
our stellar parameters. Only small changes were found. 
Temperatures changed by up to 40 K (making the stars 
warmer, although not in all cases) and the microturbu- 
lent velocity changed between 0.05 and 0.3kms~^. The 
main concern with this exercises was to have sufhcient 
number of Fel lines (between 50 and 150) and Fell lines 
(more than at least 2 to 3) available, many of which 
(especially Fell) are located between 4400 and 4650 A. 
Hence, for the near main-sequence turn-off stars cutting 
at 4750 A did not yield well-determined stellar parame- 
ters since only 20 (for G64-12) and 6 (for CD-24 17504) 
Fel and only 2 Fell (both stars) were available. 

Finally, we note that any line abundance corrections 
based on line formation under non-local thermodynamic 
equilibrium (NLTE) would increase the surface gravity. 
We made a test where we assumed a simple average cor- 
rection for Fe I of 0.2 dex, and assuming that NLTE does 



not affect Fell (e.g.. I Asplundl [20051 ) . This increases the 
gravity of about 0.5 dex. Such a large change might be 
welcome for the cooler giants when not applying our ad- 
justments, i.e. in combination with the too cool tem- 
peratures that result in too low gravities, but for near- 
main-sequence stars a large gravity increase would pose 
problems regardless. More detailed stellar parameter de- 
termination te chniques seem to be required to resolv e 
this issue (e.g. iBergemann et al.ll2012l : lLrnd et al.ll2012D . 

4. METAL ABUNDANCE RESULTS AND COMPARISON 
WITH LITERATURE VALUES 

For completeness, we present the chemical abundances 
of our sam ple stars in Ta b le [H Solar abundances are 
taken from lAsplund et al.l (J2009D . We note that while 
for the Fe abundances we employed only lines with wave- 
lengths longer than 3700 A, but for a few other elements 
(Ti, Cr, Ni) we also included a few lines at shorter wave- 
lengths for better statistics whenever they were mea- 
surable. For each star, we also compared the abun- 
dances with those of another study from the literature 
(INorris et al.l[200ll: iHoiida et al.l[200l ICavrel et al.]l200l 
lAoki et al.l 120060 . We chose those studies because they 
yielded stellar parameters for a given star that agreed 
well w i th ou r results. We took the abundances from 
iFrebell ()2010[ ) to ensure that all studies employ the same 
solar abundances. The stellar abundance differences in 
the final [X/Fe] values between our and the literature 
studies are reasonable and typically within 0.2 dex. 

Table [S] lists stellar parameters fo r an a dditional five 
stars from the sample of lCavrel et al.l ()2004[ ) for which we 
also had Magellan/MIKE spectra. We used these stars to 
further test our method and stellar parameter determina- 
tion technique. Furthermore, we add three of our sample 
stars (HD122563, BD -18 5550, CS22892-052) because 
they are also part of the Cayrel et al. sample. This allows 
up to compare our stellar param eters with the Cay rel et 
al. results as well as with those of lYong et al.l (|2013D who 
analyzed the Cayrel et al. stars as part of a sample of 
190 metal-poor stars from the literature. Figure [7] illus- 
trates the results and shows the good agreement of our 
effective temperature and surface gravity compared to 
the photometrically derived temperatures and gravities 
(a lso derived from th e FeTFell ioniz ation equilibrium) 
bv lCavrel et~all (|2004[ ) (top panel) and lYong et"all ()2013l i 
(bottom panel). 

The only star for which we determine very differ- 
ent stellar parameters than ICavrel et al.l ()2004[ ) and 
lYong et all (l2013l) is CS30325-094. We find its corrected 
temperature to be about 340 K cooler. A simple check 
of the shape of the Balmer lines in comparison with the 
stars from our sample suggests the star to be not warmer 
than 4800 K (on our adjusted scale). This remains in 
contrast with what has been found from photometry. 

A similar, albeit less severe case, is HD2796. Our cor- 
rected temperature is about 100 K cooler, but this star 
also has a somewhat higher metallicity ([Fe/H] ~ —2.5). 
If the differences between excitation temperatures and 
photometric temperatures was less severe for giants at 
higher metallicity, this could explain the discrepant tem- 
perature. However, both our initial and final stel l ar pa - 
rameter solutions as well as that of lCavrel et al.l (|2004[ ) 
do not match any position of the red giant branch on 
the isochrone. Rather, the star sits slightly above it. 
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TABLE 4 
Magellan/mike Chemical Abundances of the Sample Stars 







HD122563 




HE 1523- 


0901 






BD -18 5550 






CS22892-052 




Species 


Ige(X) 


[X/FcJ 


N 


a 


Ig^X) 


[X/h'c\ 


N 


o- 


lge(Xj 


[X/FcJ 


N 


cr 


Ig^X) 


[X/FeJ 


N 


a 


Nal 


3.95 


0.50 


2 


0.05 


3.52 


0.13 


2 


0.02 


3.45 


0.23 


2 


0.07 


3.62 


0.46 


2 


0.10 


Mgl 


5.26 


0.45 


7 


0.07 


5.08 


0.33 


6 


0.09 


4.99 


0.41 


7 


0.06 


4.85 


0.33 


5 


0.09 


Si I 


5.17 


0.45 


3 


0.06 


5.09 


0.43 


1 


0.00 


4.90 


0.41 


1 


0.00 


5.05 


0.62 


1 


0.00 


Cal 


3.87 


0.32 


25 


0.08 


3.78 


0.29 


21 


0.08 


3.72 


0.40 


25 


0.08 


3.66 


0.39 


18 


0.08 


Scl 


0.42 


0.07 


13 


0.07 


0.27 


-0.03 


12 


0.08 


0.17 


0.04 


9 


0.06 


0.03 


-0.04 


6 


0.07 


Til 


2.28 


0.12 


31 


0.07 


2.26 


0.16 


26 


0.09 


2.15 


0.22 


28 


0.09 


2.07 


0.19 


13 


0.09 


Till 


2.37 


0.21 


51 


0.10 


2.35 


0.24 


41 


0.11 


2.16 


0.23 


49 


0.11 


2.10 


0.22 


36 


0.10 


CrI 


2.56 


-0.29 


21 


0.08 


2.65 


-0.15 


14 


0.11 


2.43 


-0.19 


20 


0.03 


2.26 


-0.30 


10 


0.06 


CrII 


2.89 


0.04 


3 


0.05 


3.12 


0.32 


2 


0.12 


2.73 


0.11 


3 


0.05 


2.45 


-0.11 


1 


0.00 


Pel 


4.71 


0.00 
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0.12 


4.65 


0.00 


183 


0.12 


4.48 


0.00 


233 


0.11 


4.42 


0.00 


191 


0.12 


Fell 


4.72 


0.01 


26 


0.10 


4.65 


-0.00 


25 


0.13 


4.49 


0.01 


24 


0.08 


4.44 


0.01 


21 


0.10 


Col 


2.42 


0.22 


6 


0.10 


2.49 


0.34 


4 


0.20 


2.24 


0.27 


8 


0.10 


2.20 


0.29 


7 


0.10 


Nil 


3.56 


0.13 


24 


0.08 


3.40 


0.03 


17 


0.09 


3.35 


0.15 


22 


0.11 


3.36 


0.22 


11 


0.09 


Znl 


1.92 


0.15 


2 


0.04 


1.92 


0.20 


2 


0.04 
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0.33 


2 


0.05 


2.06 
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1 
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perhaps suggesting that the star is descending onto the 
horizontal branch. However, lYong et"al] (|2013[ ) do find 
a gravity that that matches the isochrone. Again, the 
shape of the Balmer lines suggest a cooler temperature 
than the photometrically derived temperatures. 



5. CONCLUSION 

We have presented a straight-forward method to ad- 
just spectroscopic temperatures that have been derived 
from Fe I line abundances as a function of their excita- 
tion potential. The main aim was to avoid obtaining 
too cool spectroscopic temperatures especially for upper 
red giant branch stars that are often several hundred de- 
grees cooler than photometrically derived temperatures. 
A pleasant side effect of this temperature adjustment 
is that higher surface gravities are obtained and artifi- 
cially high microturbulent velocities are avoided. Over- 
all this brings the fully spectroscopic, reddening inde- 
pendent stellar parameters in agreement with photomet- 



rically derived ones, and also with isochrones of the re- 
spective stellar metallicities. While the stellar sample is 
based on stars with —3.3 < [Fe/H] < —2.5, additional 
tests indicate that our adjustment scheme is not usable 
for stars with [Fe/H] < -4.0. 

We tested our m ethod with additional stars from the 
Cavrel et alj (I2004D sample that were also analyzed by 
Yong et al.l ()2013D and generally found agreement to 



within ~50K. The results of applying the new tempera- 
ture determination method to new samples of extremely 
metal-poor stars will be reported in a forthcoming paper. 
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Fig. 7. — Initial (small open circles), corrected (larger open 
circles) stellar parameters compared with literature values (crosses) 
of Cavrel et al. (2004) (top panel) and Yong et al.l l(20iD (bottom 
panel). Color code: green: HD2796, red: CS22873-166, black: 
CS30325-094, blue: HD186478, cyan: CS31082-001. Overplotted 
are 12 Gyr isochrones with [Fe/H] = -3.0, [Fe/H] = -2.5 and 
[Fe/H] = —2. (from left to rig ht), all having an a-enhancement of 
[a/Fe] = 0.3 HKim et al.ll200^) . See text for discussion. 
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